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We study the phenomenon of mass loss by a scalar charge — a point particle that acts a source for 
a noninteracting scalar field — in an expanding universe. The charge is placed on comoving world 
lines of two cosmological spacetimes: a de Sitter universe, and a spatially-flat, matter-dominated 
universe. In both cases, we find that the particle's rest mass is not a constant, but that it changes 
in response to the emission of monopole scalar radiation by the particle. In de Sitter spacetime, 
the particle radiates all of its mass within a finite proper time. In the matter-dominated cosmology, 
this happens only if the charge of the particle is sufficiently large; for smaller charges the particle 
first loses some of its mass, but then regains it all eventually. 

PACS numbers: 04.25.-g, 95.30.Cq, 98.80.-k 
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I. INTRODUCTION AND SUMMARY 

In this paper we study the phenomenon of mass loss 
by a scalar charge in an expanding universe. By scalar 
charge we mean a pointlike particle carrying a charge q 
that acts as a source for a massless scalar field $. The 
theory describing the particle-field system is a close ana- 
logue to standard electromagnetic theory. It comes with 
a wave equation for the field, and equations of motion 
for the particle; both equations are linear in 4>. The 
equations predict that the mass of the particle is not a 
constant, but that it changes dynamically as the particle 
moves in a curved spacetime This comes about be- 
cause the four-force produced by the scalar field, which 
is proportional to the gradient of <&, is not orthogonal 
to the particle's four-velocity. Although such an unusual 
situation was identified in the past g], it has attracted 
surprisingly little attention, and its consequences are well 
worth exploring. 

The usual notion of a particle's rest mass is that it is 
an invariant quantity that always stays constant. Fur- 
ther, one normally assumes that all electrons, say, have 
the same mass, which is then a fundamental constant of 
nature. Indeed, the whole concept of an elementary par- 
ticle tends to imply rest-mass universality, and rest-mass 
conservation. It is also usually taken for granted that 
the properties of elementary particles — such as mass, 
charge, or spin — are independent of the cosmological 
parameters. While suggestions were made to link cos- 
mology to the properties of elementary particles (notably 
the large-numbers hypothesis || 0, ||, |[ |?j ) , this idea has 
not gained much popularity. 

We show here that such a link is unavoidable for scalar 
charges — a scalar particle at rest in an expanding uni- 
verse possesses a changing mass whose dynamics is di- 



rectly coupled to the cosmology. We display this connec- 
tion for two spacetimes: de Sitter, and a spatially-flat 
matter-dominated cosmology. A de Sitter universe pro- 
vides a reasonable description of the inflationary epoch 
of our own universe, and the flat matter-dominated uni- 
verse adequately models our universe's present epoch (in 
the absence of dark energy). Thanks to the simplicity 
of these spacetimes, our mathematical treatment of the 
mass-change phenomenon is exact: our calculations in- 
volve no approximations. The phenomenon is not ex- 
clusive to the two spacetimes considered here: A simi- 
lar effect occurs in flat (1 + l)-dimensional and (2 + 1)- 
dimensional Minkowski spacetimes, and this is discussed 
in a separate paper j8|. The effect, however, does not 
occur in (3 + l)-dimensional Minkowski spacetime, nor 
in the Kerr family of spacetimes || . 

For both cosmologies we calculate the changing mass 
of a scalar charge taken to be comoving with the cosmo- 
logical fluid. We assume that the particle starts with a 
finite mass mo in the finite past, and we determine how 
the mass behaves thereafter; its behavior depends on mo 
and q, as well as the cosmological parameters. For our 
matter-dominated universe, we find that two types of be- 
haviors are possible, depending on the size of q relative 
to a combination of other parameters. In the first sce- 
nario (q small), the particle first loses a fraction of its 
mass, but it then regains it all eventually. In the second 
scenario (q large), the particle loses all of its mass within 
a finite proper time. For de Sitter spacetime, only the 
second scenario applies, and the particle loses mass at a 
constant rate. 

The mechanism behind the mass-change phenomenon 
is easy to identify: the scalar charge emits monopole 
waves, and the energy carried off by the radiation is taken 
from the particle's rest mass. (As we show in the pa- 
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per, this interpretation is precise for de Sitter spacetime, 
which admits a timelike Killing vector. It is not pre- 
cise, but still loosely correct, in the case of our matter- 
dominated universe, which does not admit a timclike 
Killing vector.) The radiative process is fueled by the 
expansion of the universe (which provides the required 
dynamics), and its monopolar nature reflects the spheri- 
cal symmetry of the problem. This implies that the mass- 
change phenomenon should not be expected for fields of 
higher spins: monopole waves could not be produced in 
such cases. Instead, the radiation would necessarily be 
associated with higher multipolc moments, and it could 
not be produced by a changing rest mass. For example, 
the four-force produced by an electromagnetic field is al- 
ways orthogonal to the four-velocity, and the rest mass 
of an electric charge is always constant. 

The classical framework adopted here does not allow 
us to predict what happens when a particle has radiated 
all of its mass. To avoid a runaway regime of negative 
rest masses, we assume that the scalar charge simply dis- 
appears when its mass drops to zero. This fix is im- 
posed without a proper justification, but it seems rea- 
sonable, and it does not appear to violate any known 
law of physics. In this regard we remark that the theory 
adopted here to describe the particle-field system does 
not automatically enforce scalar charge conservation. We 
assume, for simplicity, that q stays constant as long as 
the particle continues to exist, but that it jumps abruptly 
to zero when the particle has radiated all of its mass. It 
would be interesting to consider more sophisticated mod- 
els in which the scalar charge, as well as the mass, would 
be dynamically changing. We shall not, however, pursue 
this here. 

We note that our notion of a scalar charge — a point- 
like particle acting as a source for a massless scalar field 
- is different from other models which also describe 
scalar particles. For example, Seidel and Suen jl^] con- 
sidered a soliton which is made of a massive, sourceless 
scalar field. (For the Seidel-Suen soliton it is impor- 
tant that the field be massive: first, the frequency of 
the soliton-star oscillations vanishes without the mass 
term, and the solution becomes unstable; second, from 
the studies of gravitational critical phenomena we know 
that a massless scalar field either disperses to infinity 
or collapses to a black hole Jll[.) The Seidel-Suen soli- 
ton is very different from our notion of a scalar charge. 
It would be interesting, however, to investigate how the 
Seidel-Suen soliton behaves in a cosmological spacetime. 

We begin in Sec. II with the presentation of an action 
principle for the particle-field system, and a derivation 
of the equations of motion. In Sec. Ill we describe a 



monopole w aves. The mass loss is computed in Sec. VI 



class of cosmological spacetimes that includes the two 
cosmologies of interest to us. In Sec. IV we calculate 
the retarded Green's function for a scalar field living in 
these two spacetimes. An alternative derivation, specific 
to de Sitter spacetime, is described in Appendix A. In 
Sec. |v| we compute the field of a scalar charge at rest in 
the two spacetimes, and show that the particle radiates 



Finally, Sec. VII contains a discussion of energy conser- 
vation in de Sitter spacetime, and an examination of the 
implications of the mass-loss effect on the abundance of 
scalar charges in our own universe. We also consider 
how our purely classical treatment might be related to a 
fundamental quantum theory of radiating scalar charges. 
Throughout the paper, except when stated otherwise, we 
use geometrized units in which c = G = 1 . 



II. ACTION PRINCIPLE AND EQUATIONS OF 
MOTION 



A particle of bare mass mo and scalar charge q moves 
on a world line z a {X) in a curved spacetime with met- 
ric g a p; the world line is monotonically parameterized by 
A. The particle creates a scalar field <J>, and the dynam- 
ics of the particle-field system is described by the action 
principle 



S 



8tt 



(m - <z$) 



. . a 8Ax — z(X)) „ 



-gd^x, 



(2.1) 



where $ ;Q , = d$/dx a , z a = dz a /d\, g is the metric 
determinant, and 8^(x — z) is a four-dimensional Dirac 
distribution which satisfies J 8^{x — z)d A x — I if z a is 
within the domain of integration. We note that the action 
S is invariant under a reparameterization of the world 
line. After variation, it is convenient to set A equal to the 
particle's proper time r, which is obtained by integrating 
dr = ^ —g a pz a zP dX. We also note that the metric g a p 
does not participate in the dynamics; it is a prescribed 
tensor field in spacetime. 

Variation of the action with respect to $ produces a 
linear wave equation for the field, 



g af3 V a Vf3<S> = -47TM, 



where 



, . 8 A {x - z(t)) 
H(x) = q I 1= dr 



(2.2) 



(2.3) 



is the scalar charge density. Variation of the action with 
respect to z a (X) produces equations of the motion for the 
particle, 

™(t) — = q(g a ? + u a u^ <0 , (2.4) 

where u a = dz a /dr is the four-velocity, Du a /dr = 
the covariant acceleration, and 



m(r) = mo — q<& 



(2.5) 
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the dynamical mass of the particle. This last equation is 
equivalent to the differential statement 



-T- = -q<z>, a u a , 

dr 

which is independent of the bare mass to . 



(2.6) 



We note that Eqs. (gj)-(gj) and (JO) also follow from 
an action principle proposed by Quinn |l| : 



S' 



5i(x — z{t)) 



-m(r) g a 0u a u f} 
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dr >\/—gd x. 



(2.7) 



This action involves only the dynamical mass m(r), but 
it is not invariant under a reparameterization of the world 
line. 

The actions S and S' both produce a dynamically 
changing mass for a scalar charge. It is possible to 
construct an action that produces a consta nt m ass. 
For example, replacing (too — m Eq. (2.1) by 
too exp(— g$/mo) gives rise to the equations of mo- 
tion m Du a /dT — q(g af3 + u a vP)Q,p. The price to 
pay, however, is high, as Eq. (2.2) must now be re- 
placed by the nonlinear wave equation <7 a ^V a V = 
— 47r/zexp(— q$/too). 

have only formal valid- 



As they stand, Eqs. (2.4)-(2 



ity because the field $ derived from Eqs. (2.2) and (2.2) 
is singular on the world line. Quinn |lj] has shown that 
the singular part of the field does not affect the motion 
of the particle, which is then governed entirely by the 
smooth (or tail) part of the field. Quinn has thus pro- 
duced regularized versions of the particle's equations of 
motion: 



Du a 



1 



-q 2 (R a l3 uf l + u a R Pl u i:i u<) 
+ q 2 (g al3 + u a u f3 ) f T ~ G,p(T,r')dT' 



and 



^ = ~ q 2 R-q 2 u a £' G, a (r,r')dr'. 



(2.8) 



(2.9) 



Here, R a is the spacetime's Ricci tensor, and R is the 
Ricci scalar. The quantity G(x,x') appearing inside the 
integrals is the retarded G reen 's function associated with 
the scalar wave equation (|2.2|); it satisfies 



g a0 V a Vf3G(x, x 1 ) = -4tt 



Si(x — x 1 ) 



(2.10) 



where x is identified with z(t), the current position of 
the particle, while x' is identified with z(t'), the parti- 
cle's past position. The integrals extend over the entire 
past world line of the particle, from r' = — oo to (al- 
most) the current time, r' = t~ = r — e, where e is 



infinitesimally positive fij]. The integration is cut short 
to avoid the singular behavior of the Green's function as 
x' approaches x; it involves only the smooth part of the 
Green's function , wh ich is ofte n referred to as its "tail 
part". In Eqs. ( p^) and (Q, the four- velocity u a is 
evaluated at the current time r, and the Green's func- 
tion is differentiated with respect to x before making the 
identification x — z(t). 



We will evaluate and solve Eq. (2.9) for the dynamical 
mass of a scalar charge at rest in an expanding universe 
in this case the right-hand side of Eq. (2.8) vanishes 
and the particle follows a geodesic of the spacetime. We 
will see that solv ing Eq. ( |2.9| ) produces a regularized ver- 
sion of Eq. (2.5). 



III. COSMOLOGICAL SPACETIMES 

For simplicity we consider spatially-flat cosmologies, 
and write the metric as 



ds 2 



a 2 (r])(-dr] 2 + dx 2 + dy 2 + dz 2 ), 



(3.1) 



in terms of a conformal time r\. The spacetime is filled 
with a homogeneous fluid of density p, pressure p, and 
four-velocity 



dx a 
~d7 



= (a" 1 , 0,0,0), 



(3.2) 



where r is proper time for comoving observers, related to 
the conformal time by dr = a(rj) dr). The behavior of the 
scale factor a(r)) is governed by an energy-conservation 
equation, (a 3 p)' + p(a 3 )' = 0, in which a prime indi- 
cates differentiation with respect to to and Raychaud- 
huri's equation, a' 2 = (87r/3)pa 4 . 

We shall consider power-law cosmologies, for which the 
scale factor takes the simple form 



a(r,) = Crf, 



(3.3) 



where C and a are constants. For such cosmologies, the 
density is given by 



P = 



3a 2 
8ttC 2 rj 



1 



2a +2 ' 



the pressure by 



P 



a 



3a 



■p, 



and the Ricci scalar by 

6a(a-l) 1 



R 



C 2 



>1 



2a+2 ' 



(3.4) 



(3.5) 



(3.6) 



Two special cases, which constitute an equivalence 
class in a sense to be described below, will be of interest 
to us. The first is a matter-dominated cosmology char- 
acterized by a = 2, which produces a vanishing pres- 
sure. The second is a de Sitter cosmology characterized 
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TABLE I: The cosmological models considered in this pa- 
per. In both cases the scale factor is given by a(rj) = Crf , 
and proper time is defined by dr = a(rf) drj. For both cos- 
mologies, the table displays the parameter a, the scale fac- 
tor (in terms of both r\ and r), the relations r[rj) and ??(t), 
the density p(jf), the pressure-to-density ratio, and the Ricci 
scalar R{rj). In the spatially-flat, matter-dominated cosmol- 
ogy, both C and rj are positive, and the universe is expanding. 
In the de Sitter cosmology, both C and r\ are negative, and 
a(rj) = \C\/(—rj) also describes an expanding universe. (The 
maximally extended de Sitter spacetime also includes a pre- 
ceding contracting phase which we do not consider here.) 



IV. GREEN'S FUNCTION 



cosmology 


matter dominated 


de Sitter 


a 


2 


-1 


a(ji) 


<V 


|C|/(-»?) 


o(t) 


C(3r/C) 2/3 


\C\f>W 


r(jj) 


Cr, s /3 


-\c\H- 


T){t) 


(3r/C*) 1/3 


_ e -r/|C| 


P 


3/(2ttC*V) 


3/(8ttC* 2 ) 


p/p 





-1 


R 


12/(CV) 


12/C 2 



by a = —1, which produces a constant density and a 
pressure p = —p. In Table I we summarize the proper- 
ties of these cosmologies. 



The scalar wave equation (2.2) can be simplified if we 
introduce the auxiliary field variable ip defined by 



®(ri,x) 



1 



a{rj) 



if>(ri,x), 



(3.7) 



in which x — (x,y, z) represents the spatial coordinates. 
This new variable satisfies 



a(a — 1) 



ip = —Awa/j,, 



(3.8) 



I- d 2 /dy 2 + d 2 /dz 2 is the flat-space 
On the basis of Eq. 



(3. 



— +V 2 4 
drf 

where V 2 = d 2 /dx 2 
Laplacian operator, 
state the following interesting property: Two differing 
cosmological models with identical values for (3 = \(x{a — 
1) will produce identical reduced fields ip 1 provided that 
the source term a 3 /x is the same in both cases, and that 
ip starts with the same initial conditions. In this specific 
sense, we may say that two cosmological models with 
equal j3 are "equivalent" . For a prescribed value of (3, the 
two equivalent cosmological models will be characterized 
by a — a±, where 



a± 



1 1 rr 
= - ± -Jl 

2 2 V 



(3.9) 



Notice in particular that if (3 — 1, then a + = 2 and a_ = 
— 1, and we obtain the cosmological models summarized 
in Table I. These models are therefore equivalent in the 
sense adopted here. 



In order to find solutions to Eq. (2.1C), we factorizc 
the Green's function according to 



G(x,x') 



1 



-rg(x,x'), 



(4.1) 



which produces a reduced version of Green's equation, 

1_ V7^ J- 1 

drf- 



+ V 2 + -^)g{x, x') = -4tt(5(?7 - rj')5 3 (x - x') . 



if 



In the following two subsections we shall solve Eq. ( |4.2| ) 
using two different methods: the first relies on 
Hadamard's theory, and the second is based on a mode 
decomposition of Green's equation. An alternative 
derivation, specific to de Sitter spacetime, is described 
in Appendix A. 



A. Solution by Hadamard ansatz 

For (3 = (flat spacetime), the retarded solution to 
Eq. (O) is 



g tia \x,x') = 



S(u) 



(4.3) 



where u = rj — rj' — \ x — x' \ is retarded time; we see 
that the flat-spacetime Green's function has support only 
on the past light cone of the field point x. Relying 
on Hadamard's general theory [jl2[ we expect that 
g(x,x / ) will have support inside the light cone as well, 
and we write it as 

g(x, x') = 5 flat (x, x') + B(x, x')8(u), (4.4) 

where 9{u) is the Heaviside step function and B(x,x') a 
two-point function to be determined, but which is known 
to be s moot h when u — 0. [Because the right-hand side 
of Eq. (4.2) is nothing more than a flat-spacetime <5- 



function, there is no need to modify the S(u) part of the 
Green's function: it is simply equal to the flat-spacetime 
result. This assumption is justified by the following.] 
Substituting Eq. (Q into Eq. Q), we find that 5 flat 
takes care of the four-dimensional ^-function, and that 
the remainder vanishes as a distribution if B satisfies 



and 



(a + V)B = 



B ta {x-x') a +B- - V = O, 



(4.5) 



(4.6) 



where □ = 



-d 2 /d v 2 



is 



the flat-spacetime 
d'Alembertian op erat or and V(x) = 2/3/?y 2 . The right- 
hand side of Eq. ( f4.6| ) is a priori arbitrary, but it must 
vanish when u = 0; it is also constrained by the fact that 
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a simultaneous solution to both Eqs . (4.5) and (4.6) must 
exist. As a consequence of Eq. (4.6) and the fact that B 
is smooth at u = 0, we establish the coincidence limit 



lim B(x,x') — — V. 

x — >x' 2 



(4.7) 



This boundary condition allows us to solve Eqs. (4.5) and 
(4.6) uniquely. 

Fo r si mplicity we shall set the right-hand side of 
Eq. (4.6) to zero. As we shall see, this will produce the 
condition (3 = 1, and we shall therefore be restricted to 
the cosmological models summarized in Table I. 

Let x and x' be fixed points in spacetime, and let the 
relations x" a — x' a + (rj" — rj')n a describe a straight line 
going from x' to x as the parameter 77" ranges from 77' to 
77; the tangent vector n a — dx a /drj" is normalized so that 
n n = 1. Rewrite Eq. (4.6) in terms of the variable x" and 



note that the derivative of B in the dire ction of (x" 
is equal to 77" dB / 'drj' '. Equation (|i 
re-expressed as 



can therefore be 



V 



dB 
drf 



B 



1 



v(v"), 



(4- 



and this can be straig htfo rwardly integrated. The solu- 
tion that satisfies Eq. (4.7) is 



B(x,x') 



1 



2(77-77') 



V(ri")dri". 



(4.9) 



With V = 2/3/?7 2 we have that B(x,x') = P/(ijrj} 



Thus far we have generated a solution to Eq. (4.6 ) only. 
(Recall that we have set O = in this equation.) W e 
must now check that this is also a solution to Eq. (4.5). 



Because both V and B are proportional to (3, it is easy 
to see that this produces a constraint on the value of (3. 
As we have indicated previously, B(x,x') — (3/(r)r]') is a 
solution to both Eqs. (|4.5| ) and ( |4.6| ) if and only if j3 = 1. 

Our conclusion is that for (3 = 1 (which impli es e ither 
a = 2 or a = —1), the retarded solution to Eq. (4.2) is 



g(x,x') = 



d(rj — rj — \x — x'\) 6{r) — vj — \x 



\x — x'\ r\rj 

(4.10) 

The cosmological models summarized in Table I are 
equivalent in the sense used before, but also in the sense 
that they come with the same reduced Green's function 
g(x, x'). Because the scale factors are different, however, 
the actu al retarded Green's function G(x,x'), given by 
Eq. ( [4.l|) , takes a distinct form in each spacetime. 



B. Solution by mode-sum 

The method used in the preceding subsection to gener- 
ate the retarded Green's function was limited to the case 
(3=1. Here we describe an alternative method which 
could, if desired, be applied to any value of (3. Here we 



set (3 = \l(l + l) , and we notice that this covers cosmolog- 
ical models characterized by either a = I + 1 or a = — I. 
We will c onsider the case / = 1 in detail, and reproduce 

Eq. ra . 

We expand the reduced Green's function g(x, x') in 
terms of plane-wave solutions to Laplace's equation, 



g{x,x) = 



(2tt) 3 



~g( V ,7 1 ';k)e ik < x - x 'U 3 k, (4.11) 



and we substitute this into Eq. (4^). The result, after 
also Fourier transforming 63 (x — x'), is an ordinary dif- 
ferential equation for g(rj, rj'; k): 



1(1 + 1) 



drj 2 



g = 4n6(i] - rj), 



(4.12) 



where k = k-k. To generate the retarded Green's func- 
tion we set 



3(77, 77'; fc) = 0(77 - 77') 3(77, 77'; k), 



(4.13) 



in which we indicate that g depends only on the mo dulu s 
of the vector k. Substitution of Eq. ( |4.13| ) into Eq. fl4.12j ) 
reve als th at g must satisfy the homogeneous version of 
Eq. (gig), 



dr/ 2 



k 2 



1(1 + 1) 



.9 = 0, 



(4.14) 



together with the boundary conditions 



0(77 = 77'; k) = 0, 



< $(il = v'\k) = 4tt. (4.15) 
drj 



Substitution of Eq. (4.13) into Eq. (4.11) and integration 
over the angular variables associated with k yields 



g(x,x') 



where R 



0(v-rf) 
2n 2 R 



g(ri,r]';k)ksm(kR)dk, (4.16) 



For 1 = 0, the unique solution to Eqs. (4.14) and (4.15) 



is 



47T 

gi=o(v,V>k) = — sin(fcAr7), 



(4.17) 



where At/ = 77 — 77' . Substituting this into Eq. (|4.1( 
returns the flat-spacetime Green's function of Eq. (p 
To derive this we make use of the identity 



sin(wa:) sin(wx') duj = 5(x—x') — 5(x+x'), (4.18) 



and we note that the second <5-function is eliminated by 
the step function 9(Ar]) in Eq. ( 4.16| ). 



If I is an integer different from zero, Eq. ( 4.14 ) can 
be solved in terms of spherical Bessel functions. (If I is 
not an integer, the equation can still be solved in terms 
of ordinary Bessel functions.) It is simpler, however, to 
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generate solutions with the ladder operator L\ = —d/dr]+ 
(1 + 1) /rj. This works as foll ows: Suppose that we already 
have gi, a solution to Eq. (4.14) wit h a g iven value of I; 
then gi + i = Ligi is a solution to Eq. ( 4.14 ) with I replaced 
by I + 1. By repeat ed ap plication of the ladder operator, 
a solution to Eq. ( 4.14 ) with any integer value of I can 
be obtained from a seed solution go- 

We use this procedure to genera te a s olution gi =1 that 
satisfies the boundary conditions (4.15). For seed solu- 
tions we use the set {sin(fcA?y), cos(fcA?y)}. After applica- 
tion of L , we find that gi—\ must be given by a superpo- 
sition of the linearly independent solutions cos(fcA?7) — 
(krj)^ 1 sm(kArj) and sin(fcA?y) + (krj)^ 1 cos(kArj). The 
coefficien ts are arbitrary functions of rj , and after impos- 
ing Eqs. ( 4.15 ), we find that the appropriate combination 
is 



4-7T 



1 + ¥w) sin(fcA ^ 



Ar/ 



kr\r\ 



- cos(kArj) 



(4.19) 



Substituting this into Eq. ( [4 . 1 6| ) and using Eq. ( |4 . 1 q ) 
yields 

^(^^^ + ^^,4 (4.20) 



where 



R 



I(A V ,R) = — 



VV 



s'm(kR) 



sin(fcA?7) 



7TR J k 2 

- (kArj) cos(fcA?7)] dk. (4.21) 

We evaluate this by integrating the first term by parts; 
this gives 

2 f°° 1 

I(Ar), R) = - I - sxaikArj) cos(kR) dk 
n Jo k 

= 6(A<q-R). (4.22) 

We therefore arrive at 

S(A V -R) e(A v -R) 
gi=i(x,x) = h : , (4.23) 



R 



T)T) 



which agrees with Eq. (4.10). 

The method described here could be extended to other 
values of I, but we shall not pursue this here. 



FIELD OF A STATIONARY CHARGE 



A particular solution to Eq. (2.2) is 

$(x)= ( G(x,x') f i(x')y/^g 1 d' i x', 



(5.1) 



where? Gix.x') is the retarded Green's function of 

Eq. ( 2.1 0| ) and g' is the determinant of the metric eval- 
uated at x'\ the integration is over the entire spacetime 



manifold. In terms of the r educ ed variables ip(x) and 
g(x,x') — cf. Eqs. (5^) and (11) — we have 



ip(x) 



g(x,x')a a (r]')fx{x')(^x' 



(5.2) 



In this section we evaluate this for a stationary scalar 
charge. 

The scalar charge density of a point particle is given 
by Eq. (2.3). For a particle at rest (comoving) in an 
expanding univ erse with metri c (|3.ll ), the four- velocity 
is given by Eq. (|3.2| ), and Eq. ( [2. 3D gives a 3 (rj) fi(r] , x) = 
q8z(x), if we choose x — to represent the particle's 
position. (Because the spacetime is homogeneous, there 
is no loss of generality in this choice.) Substituting this 
into Eq. (5^2) gives 



ip(rj,x) = / qg(ri;x,r]';0)dr)', 



(5.3) 



where the reduced Green's function is given by Eq. ( 4.10 ). 

For a matter-dominated cosmology (a = 2, = 1 - 
see Table I), the integration starts at rj' — where the 



Green's function is singular, and the integral of Eq. (5.3) 
is logarithmically divergent. For a de Sitter cosmology 
(a — — 1,(3 — 1), the integration starts at r/' — —oo and 
the integral is also logarithmically divergent. To avoid 
this pathology, we assume that the scalar charge came 
into being in the finite past, and we let 



Q0(V - Vo) 



in Eq. ([5.3]), which becomes 



ip(rj, x)=q g(rj; x, r(\ 0) drjf . 



(5.4) 



(5.5) 



The replacement of Eq. (5.4) describes the sudden cre- 
ation of a scalar charge at a time r\ = rjo ^ {0, — oo}, 
and this constitutes a simple cure for the pathology of 
Eq. (^). It is also possible to let the charge adiabati- 
cally "switch on" , but this would needlessly lead to more 
complicated expressions. We recall that there is no law of 
charge conservation in this theory: a scalar charge can be 
spontaneously created provided that a sufficient amount 
energy is made avail able . 

Integratio n of E q. (5.5), with the reduced Green's func- 
tion of Eq. (4.11), is elementary, and we obtain 



ip(v, x ) = -0(r]-r-rto) 
r 



V 



r I rj — r 
— In 



(5.6) 



where r — \x\ is the coordinate distance from the origin. 
Equation ( |5.6| ) describes the reduced scalar field of a sta- 
tionary charge at x = 0; the full scalar field is given by 
$(?7, x) = jj}{r],x)/a(ri). The step function in Eq. ( |5.6| ) 
indicates that the moment of charge creation is registered 
at a time rj = tjq + r by an observer at a distance r from 
the charge: the information travels at the speed of light. 
Equation (5.6) implies that a stationary scalar charge in 
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an expanding universe radiates monopole waves: While 
the reduced scalar field ip is stationary in the immediate 
vicinity of the charge, the reorganization of the field lines 
caused by the underlying spacetime curvature (which is 
dynamical) produces radiation. The expansion of the 
universe also participates directly in the dynamics of the 
scalar field; this is reflected by the presence of the scale 
factor in the relation^ =ip/a. 

The field of Eq. (5.6) is singular at x — 0, where the 
particle is located. We define a renormalized local field 
V'rcn by first removing the singular part of ip and then 
taking the limit r — * 0. This gives 

^ren(r?) = - % - 7? ) 0(77/770). (5.7) 

We also define $> Te n(v) = ^teniv) / a (v)> the physical sig- 
nificance of this quantity will be revealed in the next 
section. 



VI. MASS LOSS 



In this section we evaluate Eq. (2.S) for a station- 
ary scalar charge in the two cosmological spacetimes de- 
scribed in Table I. For this situation, the particle's four- 
velocity is given by Eq. (3.2), and the relevant (smooth) 
part of th e Green's function is obtained from Eqs. (4.1) 
and (pLlOl): 



G 



smooth 



(x,x') 



C 2 (rpf) 



a + l ' 



(6.1) 



We recall that the scale factor is gi ven by a(rj) = Crj a , 
where C and a are constants; Eq. (6.1) is valid if and 
only if a is restricted to the values 2 (matter-dominated 
cosmology) and —1 (de Sitter cosmology). We recall also 
that for these cosmologies, the Ricci scalar is given by 
Eq. Q. 

For a stationary particle, u a G. a reduces to a~ 1 d v G, 
and the proper-time integral of Eq. ( |2.9| ) can be expressed 
as an integral over dr]' = a (rf) dr'. This yields 



dm 
dr] 



12 



(r),r)')a(r)')dT]' (6.2) 

J no 



for the rate of change of the particle's dynamical mass . 
Notice that we have incorporated the cutoff of Eq. (5.4) 
into this expression. 



A. Matter-dominated cosmology 



We now eval uate Eq. (6.2) for a — 2. Substitution of 
Eqs. (|t|) and (U) yields 



dm 
dr] 



Cr] 4 



[1-3111(77/770)], 



(6.3) 



and this can be immediately integrated: 



771(77) = m - 



Cif 



1*1(77/770), 



(6.4) 



where mo = 777(770). This result can b e re stated in terms 
of the renormalized local field of Eq. 15. 



(6.5) 



777,(77) = ™0 - ^011(77), 

which is analogous to Eq. (|2.5|) . 

In order to analyze Eq. ( |6.4| ) we express it as 



777(77) 



1 



Mx) ee f(x;c), 



777o X° 

in terms of the rescaled quantities 

x = — c= ^ 

770 ' Cm r] Q 3 ' 



(6.6) 



(6.7) 



The function f(x; c) is defined in the interval 1 < x < 00. 
It initially decreases from /(l; c) = 1 as x increases, and 
reaches its minimum value / min = 1 — c/(3e) at x = 
e i/3 ^ 1 3956. Then f(x;c) starts to increase, and it is 
eventually restored to its original value, f(x — » 00; c) — > 
1. For c > 3e ~ 8.1548, f(x; c) < in an interval around 
x = e 1 / 3 . 

^From these considerations, we conclude that for q 2 < 
3eC777o77o 3 , the particle first loses mass, but eventually 
regains all of it as 77 — > 00. For q 2 > 3eCmo77o 3 , on 
the other hand, the particle radiates all of its mass in 
a time shorter than A77 = (e 1 / 3 — 1)770 — 0.3956?7o; as 
was conjectured in Sec. I, this presumably signals the 
destruction of the scalar charge. 



B. de Sitter cosmology 

For a = —1, Gs mooth fa v') is a constant, and the in- 
tegral term of Eq. (6.2) contributes nothing to the mass 
loss. Instead, this comes entirely from the Ricci-scalar 
term, and we find 



G?777 



q 

c 2 ' 



(6.8) 



In terms of proper time t, mass is being lost at a constant 
rate. This result can also be expressed as 



m(r) = m - ^ (r - t 



(6.9) 



where too = to(to) is the initial mass of the scalar charge, 
and the proper time To is related to the conformal time 
770. (The relevant relations between proper and confor- 
mal time s are listed in Table I.) Here also we find that 
Eq. ( |6.9[ ) can be simply r esta ted in terms of the renor- 
malized local field of Eq. (5.7): 



m(r) = to - <7$ rcn (T); 



(6.10) 
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this again is analogous to Eq. (2.5). 

Our conclusion here is that the scalar charge will ra- 
diate all of its mass, at a constant rate, within a proper 
time At = C' 2 m Q /q 2 . This again signals the destruction 
of the scalar particle. 



VII. DISCUSSION: MASS LOSS IN DE SITTER 
SPACETIME 

A. Energy Conservation 

The de Sitter spacetime possesses a maximal set of ten 
linearly independent Killing vectors, and each one can 
be associated with a global conservation law. The spa- 
tial symmetries of the problem immediately imply linear 
and angular momentum conservation, but the explicit 
time dependence does not permit a hasty dismissal of 
energy conservation. Here we show that energy is glob- 
ally conserved: the energy radiated by a stationary scalar 
charge in de Sitter spacetime is exactly equal to the mass 
lost. We note that such an analysis cannot be adapted to 
the case of a matter-dominated cosmology, because this 
spacetime does not admit a timelike Killing vector. 

Let T a P be the stress-energy tensor of the particle-field 
system, and let £ Q denote the timelike Killing vector of 
de Sitter spacetime. By virtue of the equations T al3 .^ = 

and £a ; /3+£/3;a = we find that the vector j a = — T a ^ is 
divergence free. By integrating j a a — over a bounded 
four-dimensional volume V and using Gauss' theorem, 
we obtain the conservation statement 



E=- <f dS Q = 0, 

J dV 



(7.1) 



where dV is the volume's boundary, and d£ Q 
is an outward-directed surface element on dV; if 
(y 1 ,^ 2 ,?/ 3 ) are coordinates intrinsic to dV, then d£ M = 
s^ap-yidx" /dy 1 )(dx f3 /dy 2 )(dx" < /dy 3 ) d 3 y, where £ 
the totally antisymmetric Levi-Civita tensor. 

The total stress-en ergy tensor T a @ is obtained by vary- 
ing the action of Eq. (2.1) with respect to the metric. We 
express it as 



rpafi rpa(3 

1 ~ 1 field 



T 



a/3 



particle ' 



where 



is the stress-energy tensor of the scalar field, and 



(7.2) 



(7.3) 



dr (7.4) 



is the stress-energy tensor of the particle. 

We want to evaluate Eq. (fTl]) when dV consists of 
a cylindrical "tube" B surrounding the particle's world 



line, closed off by "caps" C\ and Ci (each of constant 
time) at both ends. It would be inconvenient to carry 
out this computation in the cosmological coordinates 
(rj,x,y,z), because in these coordinates the metric is 
explicitly time dependent, and the Killing symmetry is 
poorly represented. We therefore prefer to use static 
coordinates (t,r* ,9,<p), in which the metric is explicitly 
time independent. The transformation is 



r\ = — cosh(fiT*)e 



(7.5) 



(x,y,z) = sinh(KX*)e Kt (sin(9cos^>,sin#sin</>, cos (9), 

where k = 1/\C\, and in the new coordinates, the de 
Sitter metric takes the form 



*2 



ds 2 



-dt 2 + dr 
cosh 2 (ftr*) 



+ k~ 2 tanh 2 {nr* ) (d6 2 + sin 2 9d(/) 2 ). 

(7.6) 

In these coordinat es, = S a t , and the Green's function 
of Eqs. Q) and fl4.1C| ) is 



G(x,x') = KC0th{nr*)5(t-t' -r*) + n 2 6(t-t' -r*) (7.7) 

if the source point is at the spatial origin of the coordinate 
system. The field of a stationary charge at that position 
is given by 

$(t, r *) = qK 0(t - r* ~ t ) \coth(nr*) + n(t - r* - t ) , 

where to denotes the time at which the scalar charge came 
into being. Because t is proper time for an observer at 
r* = 0, the dynamical mass of the scalar particle can be 
expressed as 



m(t) = mo — q k (t — to) 



(7.9) 



where mo is the pa rticl e's mass at the time i ; this follows 
directly from Eq. ( |6.9[ ). 

We choose the closed hypersurface dV to be the union 
of a three-cylinder B described by (r* — R, t\ < t < ti), 
a spherical ball C\ described by (t = ti,0 < r* < R), 
and another ball C2 described by (t = £2 , < r* < R). 
After decomposing the total stress-energy tensor as in 
Eq. (O), we find that Eq. (O) becomes 



E — EfieldL-B] + Afield [C2] — Afield [Cl] 



-Eparticle [B] + -Eparticlc [C2] 



0, 



-^particle [Cl] 

(7.10) 



where, for example, E Rc \ d [B] = - J B T fiel |^ dT, a . The 
minus signs in front of the terms associated with C\ re- 
flect the fact that the future-directed surface element on 
Ci points within the region V, and therefore against the 
outward-directed surface element on dV . 

Omitting all calculational details, we now present our 



results for the various quantities appearing in Eq. (7.10) 
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First, we find that 



2,„2 



-Bfieid[C 2 ] = £fi e id[Ci] = -q 2 K 

r R 

x / [tanh 2 («;r*) + coth 2 (Kr*)] dr* . 
Jo 

(7.11) 

This represents the field energy enclosed within a sphere 
of coordinate radius R centered at the origin. Because the 
field configuration is singular at r* = 0, both Efi e \a[C 2 ] 
and -Efi c id[Ci] are formally infinite. A proper derivation 
of global energy conservation should therefore involve a 
regularization procedure for the field's energy. Fortu- 
nately, equality of Eneldi^] an d -Efieid[Ca] ensures that 
independently of the details of the regularization proce- 
dure, these quantities cancel each other out on the right- 
hand side of Eq. ( 7.1 0| ) . (A technical requirement is that 
the regularization procedure must not violate the time- 
translational invariance of the field energy.) Second, we 
have 



Efi c i d [B] =q Z K 2 (t 2 -t 1 ), 



(7.12) 



which represents the energy radiated by the particle dur- 
ing the time interval t\ < t < t 2 (notice that this is 
independent of R) . Third, 

-Bparticlc[C2] = m(t 2 ) , ^particle [Cl] = m(t 1 ) , (7.13) 



and -Eparticie [B] = 0: the energy carried by the particle 
is equal to the current valu e of its mas s. 

Substituting Eqs. (|7 12|) and (f7.13|) into Eq. ( |7.10| ) 
yields 



E = q 2 n 2 (t 2 - h) + m{t 2 ) - m(h) . 



(7.14) 



Using Eq. (7.9) reveals that energy is globally conserved: 
the change in the mass is exactly equal to the energy 
radiated, and E = 0. 



B. Implications for this universe 

Can the mass-loss effect be used to explain the ob- 
served absence of scalar charges? We argue here that the 
answer is in the affirmative. We will calculate an upper 
bound on the total number of scalar charges that might 
still exist today, and show that it is very small, of the or- 
der of millions per galaxy. And we will calculate a lower 
bound on the mass of a scalar charge, and show that it 
is too large to permit the production of these particles in 
today's accelerators. We assume that scalar charges al- 
ready existed at the onset of the inflationary epoch, and 
that no additional scalar charges were created since. 

The following discussion must be preceded by an im- 
portant disclaimer. The mass-loss phenomenon was in- 



vestigated in Sec. VI within a framework in which the de- 
grees of freedom associated with the charge and its scalar 



field were both described by classical physics. While we 
can hope that such a treatment may not be blatantly 
at odds with a proper quantum description, the follow- 
ing speculations regarding the fate of elementary scalar 
charges must be regarded as more conjecture than defi- 
nite prediction. In this regard we can be encouraged by 
the fact that the classical treatment of electromagnetic 
radiation reaction is often not a bad approximation to 
the correct quantum description |14| ]. 

A fundamental quantum theory of radiating scalar 
charges might be based on the following considerations. 
The first ingredient would be a complex scalar field 
^particle whose quantum excitations would give rise to 
massive particles that carry a scalar charge. The mass 
of the particles would be identified by locating the poles 
of the Feynman propagator in momentum space |l5j . In 
the absence of radiative corrections, and in flat space- 
time, these particles would have a mass that would stay 
constant, and this mass would directly correspond to the 
mass parameter of the field's Lagrangian. The physics of 
these elementary excitations might, however, be substan- 
tially different in the curved spacetime of an expanding 
universe. For example, Redmount ]l6| has shown that 
in the context of a real scalar field in de Sitter space- 
time, the particle energies are not constant — they oscil- 
late at late and early times — and do not correspond to 
the Lagrangian's mass parameter. Although Redmount 's 
particles do not carry a scalar charge, his results clearly 
suggest that the quantum physics of a scalar charge in an 
expanding universe, even in the absence of radiative cor- 
rections, might be considerably richer than the classical 
description provided in this paper. 

The second ingredient involved in a fundamental the- 
ory of radiating scalar charges would be another scalar 
field ^radiation coupled to the first to allow the scalar par- 
ticles to radiate. This new field would be real and mass- 
less, and it would generate the radiative corrections that 
have so far been missing in our description. These would 
modify the energies of the scalar charges with respect to 
the (already complicated) free-field behavior. While this 
fully quantum description of a radiating scalar charge 
would undoubtedly be richer and more interesting than 
the classical treatment provided in this paper, we can 
hope that our classical considerations will not lead us 
too far astray — they should indeed provide us with a 
useful approximation. In this spirit we shall pursue our 
speculations regarding the fate of scalar charges in our 
own inflationary universe. 

For concreteness we assume that the time scale asso- 
ciated with the inflationary epoch is t c = \C\ ~ 10~ 34 s. 
This produces a distance scale r c = ct c ~ 3 x 10~ 24 cm, 
a density scale p c = 3/ (8irGt c 2 ) ~ 2 x 10 74 g/cm 3 , and a 
mass scale m c = 47rp c r c 3 /3 ~2x 10 4 g. For the purpose 
of this discussion we reintroduce the speed of light c and 
the gravitational constant G, which were both previously 
set equal to one. 

We assume that a number N of scalar charges, all of 
the same mass toqi are created prior to the onset of infla- 
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tion, and th at d uring inflation, their mass decays accord- 
ing to Eq. (7.S). An upper bound on N is obtained by 
noting that at the onset of inflation, iVm-o cannot exceed 
m c , the total mass contained in (what will become) the 
observable universe. Thus, N < m c /rriQ. 

Whether or not the scalar charges will survive the in- 
flationary epoch depends on the relation between mo and 
its duration At: if the particles are created heavy, and if 
inflation doesn't persist for too long, then every particle 
present initially will still be present after inflation. With 
our objective to produce an upper bound for N, we as- 
sume that the scalar charges do survive the inflationary 
epoch, and we use Eq. (7.9) to derive a lower bound on 
the initial mass Too- After inserting the appropriate fac- 
tors of c and normalizing the (unknown) scalar charge q 
to the electron's charge e, we obtain 



mo > 



At 
17 



(7.15) 



For an inflationary epoch that persists for approximately 
60 e-folding times, we have that m > 5 x 10~ 15 (g/e) 2 g. 
If q/e is of order unity, this mass is larger by 6 orders of 
magnitude than the energy currently available at particle 
accelerators. 

Combining this lower bound on mo with our previous 
result for N, we find that N < 4 x 10 18 (e/g) 2 , so that 
the total number of scalar charges present today cannot 
exceed 10 7 (e/g) 2 per galaxy. So we see that unless e/q 
is huge, mo is naturally extremely large, N is naturally 
extremely small, and the prospect of observing a scalar 
charge today is extremely limited. 



Acknowledgments 

Our result for the retarded Green's function of a scalar 
field in de Sitter spacetime, and the realization that scalar 
charges in de Sitter spacetime must have a variable mass, 
were arrived at independently by Amos Ori [ [L7| . We 
thank him for discussions, as well as Karel Kuchaf, Lee 
Lindblom, and Kip Thorne. L.M.B. wishes to thank the 
Technion Institute of Theoretical Physics for hospitality. 

E.P. was supported by the Natural Sciences and En- 
gineering Research Council (Canada). A.I.H. was sup- 
ported by a SURF fellowship. L.M.B. was supported by 
NSF grants PHY-0099568 (Caltech) and PHY-9734871 
(University of Utah) . 



APPENDIX A: ALTERNATIVE DERIVATION OF 
THE SCALAR FIELD GREEN'S FUNCTION FOR 
DE SITTER SPACETIME 

In this Appendix we present an alternative derivation 
for the retarded Green's function of a scalar field in de 
Sitter spacetime. 

First, we observe that because of the four-dimensional 
homogeneity of de Sitter spacetime, the Green's function 



can only be a function of the invariant distance between 
two points. Next, we adopt coordinates that make it easy 
for us to use this observation. Specifically, we write the 
de Sitter metric in the form 

ds 2 = -dT 2 + n- 2 sinh 2 (nT) ( dr 2 + sinh 2 r dtt 2 ) , (Al) 
where k = | C | 1 and dtt 2 is the standard metric on the 



unit two-sphere. In the coordinates of Eq. (Al), the spa- 



tial sections of de Sitter spacetime are open hyperboloids 
of constant (negative) curvature. This form of the metric 
can be obtained from the embedding relations 



Xq = k~ x sinh(Kr) coshr 

X\ = k" 1 cosIi(kt) 

X2 = K~ sinh(/€T) sinh v cos 9 

x% = ft -1 sinh(Kr) sinh r sin 9 cos < 

X4 = sinh(Kr) sinh r sin 9 sin <j 



(A2) 



which describe de Sitter spacetime as the hypersurface 

" 2 in a five-dimensional flat 



■x 2 



■ x 2 



■ x\ 



spacetime. The coordinate ranges are — 00 < r < 00, 
< r < 00, < 9 < 7T, and < (j) < 2tt. Notice 
that these coordinates do not cover the entire de Sitter 
manifold. They go bad onr = and r = 00, which corre- 
spond to two null rays originating from the origin. More 
importantly, r = is a coordinate singularity: all the 
points (with finite coordinate values) (r = 0, r, 9, <ft) are 
the same physical point on the manifold ]l7| ]. This can 
most easily seen by direct substitution in the coordinate 
transformation of Eq. (A2). 

These coordinates are particularly convenient, because 
surfaces of constant geodesic distance from the origin are 
also surfaces of constant r. We put the source at the ori- 
gin of the coordinates without loss of generality (because 
of the homogeneity of de Sitter spacetime). Next, we 
will show that the proper time along a geodesic between 
the points Pq at (0, 0, 0, 0) and Pi at (r, r, 9, <p) is simply 
t for all finite values of (r,9,4>). This is most elegantly 
accomplished by a method suggested to us by Ori p7| . 
From the previous result that all the points with r = 
(and with finite coordinate values) are the same point, Pq 
may be represented by the coordinates (0,r,6,4>), from 
which a constant (r, 9, 4>) curve (which is a geodesic) may 
be extended to the point Pi at (r,r,6,(f>). Next, observe 
that the proper time along geodesies of constant (r, 9, (j>) 
is simply the coordinate time. It then follows that the 
geodesic distance between Pq and Pi is r. The retarded 
Green's function sourced at (0,0,0,0) is thus a function 
only of the coordinate time for all evaluation points, that 
isG(x,0) = G(T,0) = G(T). 

Next, using the previous result that G = G(r), we 
use Green's equation to obtain an ordinary differential 
equation for G(t). Specifically, we find that 



sinlT 3 (KT)<9 T [sinh 3 (Kr)5 r G(T)] = 4tt 



<5 4 (x) 



(A3) 



Inside the source's light cone, the Green's function is a 
linear combination of two linearly-independent solutions 
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for the corresponding homogeneous equation. These two 
solutions are just 



Gi(r) = 1 and G 2 (r) = 



cosIl(kt) 1 cosh(ftr) — 1 
+ 2 n cosh(/tr) + 1 



sinh (kt) 
_ (A4) 
Hadamard's theory Q requires that G(r) = ciG\{t) + 
c 2 G 2 (r) be an analytic function. However, G 2 is not 
analytic: approaching r = it behaves like G%{t) = 
1/(2k 2 t 2 ) + ln(Y)/2 + 0(1). Hence, c 2 = 0, and we con- 
clude that the Green's function is a constant inside the 
source's light cone. 



The coordinates (Al) are inconvenient for finding the 
value of that constant, because of their bad behavior on 
the light cone. Because the Green's function is a scalar, 
our conclusion that it is a constant inside the light cone is 
unchanged when we switch to other coordinates. As was 



discussed in Sec. IV, it may be directly observed from 



the form of Green's equation in conformal coordinates 
that the support on the source's light cone is identical to 
the corresponding flat-space sup port. Transforming to 
the static coordinates of Eq. (7.6), 



G(x, 0) = KC0th( K r*)6(t - r*) + K6(t - r*), (A5) 

where K is yet to be determined. 

Substituting this into the homogeneous Green's equa- 
tion in these coordinates, 



[cosh 2 (/ 



2KC0th(Kr*)<9 r ..J G(t,r*) = 0. 

_ (A6) 

fixes K = k 2 . With this identification, E q. (Aq ) is equiv- 
alent to the Green's function of Eqs. (4.1) and (4.1C) 
adapted to de Sitter spacetime. 
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